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1  Introduction 


Most  investigations  of  internal  phosphorus  (P)  loading  in  lakes  have  con¬ 
centrated  on  P  rdease  from  profundal  sediments  under  anaerobic  xonditions 
(e.g.  Riley  and  Prepas  1984,  Ntimberg  1987).  However,  recent  studies  have 
shown  diat  P  can  also  be  released  from  littoral  sediments  under  aerobic  condi¬ 
tions  (Twinch  and  Peters  1984,  Drake  and  Heaney  1987).  Phosphorus 
rdeased  during  the  senescence  of  aquatic  macrophytes  constitutes  another 
potential  source  of  internal  P  loading  from  the  littoral  zone  (Barko  and  Smart 
1980,  Carpenter  1980).  Since  internal  P  loading  from  littoral  macrophyte 
beds  may  significandy  influence  the  P  economy  of  aquatic  systems,  particu¬ 
larly  small  lakes  (Prentki  et  al.  1979),  it  is  important  to  evaluate  horizontal  P 
exchanges  b^een  littoral  and  pelagic  zones  of  these  systems  (James  and 
Barko  1991). 


We  rqmrt  the  results  of  seasonal  investigations  of  P  exchange  between  the 
littoral  and  pelagic  zones  of  Eau  Galle  Reservoir,  Wisconsin,  driven  by  night¬ 
time  convective  circulation.  Because  of  the  importance  of  P  gradients  in 
affecting  P  flux  rates  during  these  periods  of  circulation,  we  examine  profiles 
of  P  both  in  the  sediment  and  in  the  overlying  water  column.  Rates  of  P 
release  from  littoral  sediments,  affected  by  both  pH  and  oxygen  as  a  function 
of  macrophyte  metabolism,  are  estimated  from  results  of  laboratory  incubation 
studies  employing  the  use  of  intact  sediment  cores.  The  objective  of  this 
rqwrt  is  to  expand  earlier  findings  (James  and  Barko  1991)  with  greater  atten¬ 
tion  to  both  P  release  from  littoral  sediments  and  climatic  and  limnological 
conditions  influencing  littoraii)elagic  P  dynamics  during  nighttime  convective 
circulation. 
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2  Methods 


Phosphorus  Release  from  Littoral  Sediments 

Rates  of  soluble  reactive  P  (SRP)  release  from  littoral  sediments  were 
measured  in  the  laboratory  at  pH  values  ranging  from  8.0  to  10.0  under  both 
a«:obic  and  anaerobic  conditions,  using  intact  sediment  cores  collected  from 
macrophyte  beds  within  the  littoral  zone.  Three  rqilicate  cores,  each  contain¬ 
ing  a  sediment  volume  of  330  cm^  (6.S  cm  in  diameter,  10  cm  high),  were 
collected  in  acrylic  sediment  core  liners.  The  overlying  water  was  carehilly 
drained  from  each  sediment  core  and  filtered  through  a  glass-fiber  filter 
(Gelman  A-E),  with  300  ml  then  siphoned  back  onto  the  sediment  without 
causing  sediment  resuspension. 

The  sediment  incubation  systems  consisted  of  sediment  and  overlying  water 
contained  in  core  liners  that  were  sealed  with  rubber  stoppers.  These  incuba¬ 
tion  systems  were  placed  in  a  darkened  environmental  chamber  (Psychro- 
therm)  and  incubated  at  20  ”C  for  1  to  2  weeks. 

The  oxidation-reduction  environment  in  eadi  system  was  controlled  by 
gently  bubbling  eidier  air  (aerobic)  or  nitrogen  (anaerobic)  through  an  air 
stone  placed  just  above  the  sediment  surface.  Bubbling  action  ensured  com¬ 
plete  mixing  of  the  water  colunm  but  did  not  disrupt  the  sediment.  An  auto¬ 
mated,  recording  pH  controller  (New  Brunswick)  was  used  to  maintain  the  pH 
at  a  constant  value  in  each  incubation  system  through  additions  of  either  sul¬ 
furic  acid  (0.1  M)  or  sodium  hydroxide  (0.1  M). 

Water  samples  for  SRP  determinations  were  collected  daily  from  the  center 
of  each  incubation  system  using  an  acid-washed  syringe  and  were  immediately 
filtered  dirougfa  a  0.4S-/tm  membrane  syringe  filter.  The  water  volume 
removed  from  each  system  during  sampling  was  replaced  by  addition  of  dis¬ 
tilled  water  adjusted  to  the  oxidation-r^uction  and  pH  conditions  of  the  incu¬ 
bation  system.  These  volumes  were  accurately  measured  to  adjust  SRP  for 
dUution  effects.  Rates  of  SRP  release  from  littoral  sediments  (mg  m'^  day'*) 
were  calculated  as  the  change  in  SRP  in  the  overlying  water  divided  by  time 
and  the  area  of  the  sediment  core  liner. 
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P  Concentrations  in  Interstitiai  Water  of  Littoral 
Sediments 


Three  replicate  sediment  peckers,  similar  to  those  described  in  Shaw  and 
Prq>as  (1989),  were  placed  in  the  sediments  at  station  2  (Figure  1)  in  May, 
July,  and  August  19^  to  examine  seasonally  the  vertical  patterns  of  P  concen¬ 
tration  in  the  interstitial  water  of  littoral  sediments  in  macrophyte  beds.  Acry¬ 
lic  sediment  peq)ers  contained  six  chambers  spaced  at  2-cm  intervals. 


Figure  1 .  Morphometric  contours  Im)  and  station  locations  in  the  north¬ 
west  bay  region  of  Eau  Galle  Reservoir.  Shaded  area  represents  regions  of 
macrophyte  growth  (littoral  zone) 
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Chambers  were  filled  with  distilled  water  purged  with  nitrogen  to  remove 
dissolved  oxygen  (DO).  Nucleopore  dialysis  membranes  (1.0-ftm  pore  size) 
were  secured  onto  the  water-filled  chambers  with  an  acrylic  plate.  Immedi¬ 
ately  thereafter,  the  peepers  were  placed  in  a  nitrogen-purged  distilled  water 
bath  to  maintain  anoxic  conditions  during  peeper  deployment.  The  peepers 
were  driven  6  cm  into  the  sediment,  which  allowed  exposure  of  four  chambers 
to  the  interstitial  water  of  the  littoral-sediments  and  two  chambers  to  the  water 
column  immediately  above  the  littoral  sediment  interface.  Peepers  were  equil¬ 
ibrated  for  14  days  (Shaw  and  Prepas  1989)  before  collection. 

Upon  retrieval  of  the  peepers,  water  samples  from  each  chamber  were 
collected  with  acid-wash^  SO-ml  syringes,  ^en  filtered  inunediately  through  a 
0.4S-Mnt  membrane  filter  and  sealed  in  a  small  acid-washed  vial  with  minimal 
exposure  to  air.  The  sarnfftes  were  analyzed  for  SRP  within  2  hr  of  retrieval. 

Fickean  difftisional  fluxes  (J)  of  P  across  the  littoral-sediment  interface 
were  calculated  (from  Berner  1980),  as 

j  =  e-j  .  dcidx 


where 

♦  =  porosity  (0.8)  of  the  sediment 

D  =  areal  s^iment  diffusion  coefficient  (1.24  x  lO"*)  (Quigley  and 
Robbins  1986) 

0‘^  =  tortuosity  of  sediment  (f"®  *)  (Berner  1980) 
dCldx  =  P  gradient  across  the  littoral-sediment  interface 

Phosphorus  gradients  across  the  littoral-sediment  interface  were  calculated 
from  P  concentrations  inunediately  above  and  below  (2-cm  interval)  the 
littoral-sediment  interface. 


Limnological  Profiling 


Six  stations  were  established  along  a  transect  located  in  the  northwest  bay 
region  of  the  reservoir  (Figure  I)  for  monitoring  seasonal  changes  in  limno¬ 
logical  conditions  pertinent  to  the  investigation.  Stations  1  and  2  (0.6S  m 
deep)  were  located  within  the  macrophyte-occupied  littoral  zone;  stations  3 
and  4  (1.0  and  2.0  m  deep,  respectively)  were  located  near  the  littoral-pelagic 
interface;  and  stations  S  and  6  (4.0  and  7.0  m  deep,  respectively)  were  located 
within  the  pelagic  zone.  Vertical  profiles  of  pH,  DO,  total  P  (TP),  and  SRP 
were  obtained  at  these  stations  over  biweekly  intervals  from  June  through 
early  September  1989.  All  data  reported  here  derive  from  samples  taken  and 
measurements  made  between  2000  and  2200  hr,  a  period  occurring  typically 
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at  the  onset  of  nighttime  water  column  cooling  during  the  summer  in  Eau 
Galle  Reservoir. 

Dissolved  oxygen  and  pH  were  measured  with  a  Hydrolab  Surveyor  H, 
which  was  precalibrated  with  buffers  and  Winkler  DO  determinations  (Ameri¬ 
can  Public  Health  Association  1985).  In  situ  profiles  of  pH  and  DO  wore 
detomined  at  25-cm  d^th  intervals  at  stations  1-3  and  at  50-cm  dq)th  inter¬ 
vals  at  stations  4-6  from  the  reservoir  surface  to  within  IS  cm  of  the  sediment 
surface. 

Samples  for  both  TP  and  SRP  analyses  were  collected  with  pnoimatically 
driven  close-interval  syringe  samplers,  as  described  in  James  and  Barko 
(1991).  Filtration  for  SRP  was  accomplished  in  situ  by  attaching  0.4S-/im 
membrane  syringe  filtm  to  syringes.'  At  the  littoral  zone  stations  (1-3),  sam¬ 
ples  wo'e  collected  over  S-cm  intervals  between  5  and  20  cm  above  die  sedi- 
m^t,  and  at  12'  to  2S-cm  intervals  diereafter  to  the  reservoir  surface.  At  the 
pelagic  zone  stations  (4-6),  samples  were  collected  over  2S-cm  intervals  from 
5  cm  above  the  sediment  to  the  reservoir  surface.  TP  was  analyzed  on  a 
Technicon  Auto- Analyzer  11  following  digestion  with  potassium  persulfate 
(American  Public  Health  Association  1985).  SRP  was  analyzed  colorimetri- 
cally  on  a  Pwkin-Elmer  Lambda  3b  spectrophommeter  using  the  ascorbic  acid 
mediod  (American  Public  Health  Association  1985).  The  detection  limit  for 
all  TP  a^  SRP  analyses  was  0.005  mg  L  '. 

Exchanges  Between  Zones  During  Nighttime 
Cooling 

Water  temperature  was  measured  every  30  min  at  stations  1,  2,  5,  and  6 
(Figure  1)  from  May  through  September  1989  using  recording  thermistor 
strings  (OmniData  International)  calibrated  to  die  nearest  0.1  °C.  These  mea¬ 
surements  were  made  from  the  surface  dov.'award  at  0.25-m  dqith  intervals  at 
stations  1  and  2,  at  0.5'm  depth  intervals  at  station  5,  and  at  1.0-m  depth 
intervals  at  station  6.  Wind  speed  (r^resmting  an  average  of  5-min  intervals 
over  30  min)  and  air  temperature  (OmniData  International)  were  recorded  at 
station  3  every  30  min. 

Water  temperatures  were  integrated  vertically  over  the  entire  water  column 
of  the  littoral  zone  (stations  1-2)  and  over  the  upper  1-m  water  column  of  the 
pelagic  zone  (stations  5-6)  to  examine  the  development  of  horizontal  gradients 
in  water  temperature  between  the  littoral  and  pelagic  zones  during  nighttime 
convective  cooling.  Horizontal  water  temperature  gradients  were  calculated  as 
the  difference  between  dqith-integrated  water  temperatures  in  the  littoral  and 
pelagic  zones  for  each  0.5-hr  interval.  Mean  negative  horizontal  temperature 
gradients  w^e  calculated  for  periods  in  which  d^tfa-integrated  water  tempera¬ 
tures  in  die  littoral  zone  wm-e  cooler  than  those  in  the  pelagic  zone. 
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Hydraulic  exchange  rates  during  nighttime  cooling  were  determined  from 
patterns  of  dye  dispersion  between  the  littoral  and  pelagic  zones.  A  transect 
was  established  in  the  vegetated  region  of  the  northwest  bay,  along  the  slope 
of  the  basin,  using  posts  driven  into  the  sediment  at  10-m  intervals  (Figure  1). 
Buoys  were  placed  in  the  pelagic  zone  at  10-  to  20-m  intervals  to  extend  the 
transect  into  the  open  water.  All  stations  and  distances  along  the  transect 
were  surveyed  using  a  transit. 

Hydraulic  exchange  rates  were  determined  at  fixed  biweekly  intervals 
under  prevailing  climatic  conditions  from  June  through  August  1989  using 
Rhodamine  WT  red  fluorescent  dye  (Crompton  and  Knowles,  Inc.)  as  a 
tracer.  Methods  for  determining  hydraulic  exchange  rates  from  dye  dispersion 
were  similar  to  those  described  in  James  and  Barko  (1991).  The  dye  was 
injected  homogeneously  in  the  littoral  waf^r  column  (1.0  m  deep)  near  the 
outer  edge  of  the  macrophyte  beds  interface  (Figure  1)  at  2100  hr  during  each 
monitoring  period.  Both  ^e  forward  (into  the  pelagic  zone)  and  reverse  (into 
the  littoral  zone/macrophyte  beds)  edges  of  the  dye  cloud  were  tracked 
between  0600  and  1200  hr  the  next  morning  with  a  Turner  Designs  fluorom- 
eter  (10-005R),  calibrated  with  dye  standards.  Hourly  volumetric  flow  rates 
(m^  hr‘)  into  the  littoral  and  pelagic  zones  were  calculated  as  the  product  of 
flow  velocity  (distance  of  the  edges  at  timej  minus  the  distance  at  time, 
divided  by  timej-time,),  vertical  expanse  of  the  dye  cloud  edges,  and  horizon¬ 
tal  length  of  the  littoral-pelagic  interface  (about  1,(X)0  m)  divided  by  time 
(James  and  Barko  1991). 

TP  exchange  rates  (mg  hr*)  into  the  littoral  or  pelagic  zone  were  cal¬ 
culated  as  the  product  of  hourly  volumetric  flow  rates  and  depth-integrated  TP 
in  the  region  of  flow  divided  by  the  reservoir  surface  area.  TP  from  the 
bottom  0.25  m  of  the  littoral  zone  was  assumed  to  move  with  dye  into  the 
pelagic  zone,  while  TP  from  the  top  0.75  m  of  the  pelagic  zone  was  assumed 
to  move  with  dye  into  the  littoral  zone.  These  assumptions  are  based  on 
patterns  of  dye  dispersion  observed  earlier  in  this  reservoir  (James  and  Barko 
1991),  as  well  as  in  the  current  inv^tigation.  Net  TP  flux  rates  (mg  m'^  hr  *) 
were  calculated  as  the  difference  between  TP  exchange  rates  into  littoral  and 
pelagic  zones. 


6 


Chapter  2  Methods 


3  Results 


Littoral'Pelayic  P  Gradients 

Rates  of  SRP  release  from  littoral  sediments,  measured  in  the  laboratory  at 
20  °C,  increased  with  increasing  pH  under  both  aerobic  and  anaerobic  condi¬ 
tions  figure  2).  Under  anaerobic  conditions  at  lower  values  of  pH,  sedi¬ 
ments  exhibited  hi^er  rates  of  SRP  release  dian  under  aerobic  conditions. 
Howevn,  near  pH  10,  sediments  exhibited  similar  rates  of  SRP  release  under 
aerobic  and  anaerobic  conditions. 

Ranges  in  rates  of  SRP  release  from  littoral  sediments  in  die  reservoir  were 
estimated  during  die  poiod  June  through  August  using  the  linear  relationships 
establi^ed  between  pH  and  rates  of  SRP  rdease  from  littoral  sedimoits  mea¬ 
sured  in  the  laboratory  (Figure  2).  Predominandy  aerobic  conditions  (DO  > 
1.0  mg  L*')  prevailed  during  the  evening  hours  tlvoughout  the  summer  in  the 
littoral  zone;  however,  anoxia  may  have  developed  direcdy  at  the  sediment 
surface,  particularly  during  the  late  ni^t  and  morning  hours.  Mean  daily 
temperatures  above  the  littoral  sediments  wer.$  20  "C  or  greater  (mean  =  21.S 
±  0.03  standard  error,  SE)  on  85  percent  of  die  days  measured  during  this 
period.  Mean  pH  in  bottom  waters  of  the  littoral  zone  ranged  from  >9.0 
(range  =  9.1  to  9.4)  in  June  and  July  to  8.8  in  August.  Estimated  rates  of 
SRP  release  from  littoral  sediments  ranged  from  1.9  to  5.0  mg  m'^  day  ' 

(mean  =  3.6  mg  m'^  day'  ±  0.5  SE)  undo^  aerobic  conditions  and  from  5.5 
to  7.6  mg  m'^  day'  (mean  =  6.7  mg  m'^  day'  ±  0.3  SE)  under  anaerobic 
condittons. 

SRP  in  the  interstitial  water  of  littoral  sedimoits  increased  significandy  at 
most  dqiths  below  the  sediment  interface  between  late  May  and  August  (Fig¬ 
ure  3).  Significant  differences  in  SRP  in  the  interstitial  water  as  a  function  of 
sediment  dqptfa  were  not  observed,  however,  due  to  high  variability  among 
replicates  (analysis  of  variance,  SAS  1985).  SRP  gradimts  across  the  littoral- 
sediment  inter^ice  ranged  between  3.4  and  7.3  x  10*  mg  m**  (mean  =  5.0 
X  10*  mg  m**  ±5,1(X)  SE),  resulting  in  a  mean  diffusional  flux  of  3.6  mg  m'^ 
day'  (±0.4  SE)  for  midsummer.  This  value  was  in  ^eement  with  the  mean 
rate  of  SRP  release  from  littoral  sedimoits  estimated  for  aerobic  conditions 
under  high  ambioit  pH. 
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Rgure  2.  Variations  in  rates  of  SRP  release  from  littoral  sediments  as  a 
function  of  pH  under  aerobic  and  anaerobic  conditions  for  undisturbed  iabo* 
ratory  sediment  incubation  systems 

Marked  horizontal  and  vertical  gradients  in  TP  developed  in  the  littoral 
zone  from  May  through  August  (Figure  4),  and  SRP  accounted  for  4S  to 
90  percent  of  the  measured  TP.  TP  was  generally  mudi  lower  in  the  q)ilim- 
nion  of  the  pelagic  zone  dian  in  the  littoral  zone.  Elevated  TP  at  metalimnetic 
depths  (2  to  4  m)  in  die  pelagic  zone  in  early  July  and  early  Sqitember  coin¬ 
cided  with  intrusion  of  storm  inflows  from  Ae  Ew  Galle  River  (James  and 
Barko,  unpublished  data).  TP  also  increased  in  the  hypolimnion  of  die  pelagic 
zone  from  June  through  August.  TP  gradients  both  in  the  littoral  and  pdagic 
zones  diminished  in  Sqitember  as  water  tempo’atures  decreased  and  the  mac¬ 
rophyte  beds  entered  soiescence. 


Thermal  Gradients  and  Convective  Circulation 

Mean  negative  horizontal  temperature  gradients  Gittoral  zone  wato-  tem¬ 
perature  <  pelagic  zone  water  temperature)  generally  increased  during 
declines  in  mean  daily  air  temperature  (Figure  5).  However,  this  reladon^ip 
was  not  statistically  significant.  Mean  negative  horizontal  ttmipo^re  gradi¬ 
ents  develt^  on  72  pwcent  of  the  midsummer  days  (June-August,  Figure  S) 
and  averag^  0.41  **C  (±0.04  SE).  During  seasonal  warming  in  late  May  and 
June,  mean  n^ative  horizontal  t«nperature  gradients  developed  only  during 
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Figure  3.  Seasonal  variations  in  mean  1  SE;  n  =  3)  SRP  concentrations 
in  the  interstitial  water  of  littoral  sediments  and  overiying  water  column 
measured  at  station  2.  Different  letters  indicate  significant  differences 
(p  <  0.05)  between  dates,  based  on  Duncan's  multiple  range  analysis  (SAS 
1985) 

periods  vriien  cold  fronts  passed  diiough  die  region.  Thus,  only  31  percent  of 
the  days  monitored  in  May  and  43  percent  of  those  monitored  in  June  exhib¬ 
ited  mean  negative  horizontal  temperature  gradients.  In  July  and  August, 
however,  mean  n^ative  horizontal  temperature  gradients  developed  on  81  and 
95  pwcent  of  the  days,  respectively. 

During  each  dye-tracking  investigation,  the  littoral  zone  exhibited  cooler 
dqith-integrated  water  tenperatures  than  the  pelagic  epilimnion  during  the 
ni^t  and  morning  hours  (Figure  6).  Mean  negative  horizontal  tempwature 
gradients  during  these  hours  were  greatest  in  early  July  and  early  August 
(Table  1).  During  midday  hours  of  each  investigation,  the  opposite  pattern 
usually  occurred,  with  dqith-integrated  water  tenperatures  in  Ae  littoral  zone 
exceeding  those  in  the  pelagic  zone  (Figure  6).  Wind  speed  was  goierally 
greatest  during  die  day,  and  declined  to  <2  m  sec*’  at  mght  during  eadi 
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Figure  4.  Seasonal,  vertical,  and  horizontal  variations  in  TP  concentration  in  the  littoral  and 
pelagic  zones 
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Figure  5.  Seasonal  variations  in  mean  daily  air  temperature  (a)  and  mean 
negative  horizontal  temperature  gradient  (b) 

investigation.  Wind  speed  was  near  zero  at  night  in  late  June,  late  July,  and 
early  August  (Figure  6). 

Dye  dispersed  as  a  bottom  current  from  the  littoral  zone  to  the  pdagic  zone 
at  ni^t  (Figure  7).  Intrusion  of  littoral  water  into  the  pelagic  zone  occurred 
as  an  interflow,  consistently  confined  to  0.2S  m  in  vertical  expanse.  The 
tenq>erature  of  this  interflow  coincided  with  temperatures  measured  in  the 
littoral  zone  (Figure  7).  Due  to  seasonal  variations  in  pelagic  thermal  stratifi¬ 
cation,  the  dqpth  of  interflow  decreased  to  a  minimum  in  July,  then  increased 
in  August. 

Hourly  volumetric  flow  rates  into  both  die  littoral  aiKl  pelagic  zones  were 
similar  throughout  the  study,  indicating  an  approximate  balance  in  water 
exchange  between  the  two  zones  (Table  1).  Hourly  volumetric  flow  rates 
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Figure  6.  Diet  variations  in  air  temperature,  depth-integrated  littoral  and  pelagic  water  tem¬ 
perature,  and  wind  speed  (near  the  bottom  of  each  panel)  on  1 4  June,  28  June,  1 3  July, 
27  July,  9  August,  and  23  August  1 989 
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Table  1 

Eatimation  of  Houriy  Volumetric  Flow  Rates  (Calculated  from  Dye 
Movement  Patterns)  and  Mean  Negative  Horizontal  Temperature 
GracDents  During  Periods  of  Nighttime  Convective  Circdation 


Oat* 

Diraction 

of 

Mowaraant 

O 

III 

Voiumatrie 
Flow  Rat* 
m*iir’ 

X  10* 

Mean 

Vciumatric 
Flow  Rat* 
m*hr-’ 

Maan  l±1  8E) 
Hcriiotitai 
Tamparatur* 
GraJiant’ 

•C 

13-14  Junes 

Palagtc 

Uttoral 

0.0017 

0.0006 

1.6 

1.4 

1.4 

-0.41*  (0.061 

27-28  Jun  89 

Magic 

0.0012 

1.1 

Littoral 

0.0004 

1.1 

1.1 

-0.18*  (0.02) 

12-13  Jui  89 

Magic 

0.0038 

3.2 

Uttoral 

0.0012 

3.1 

3.2 

-0.68*  (0.071 

28-27  Jul  89 

Magic 

0.0010 

0.9 

Littoral 

0.0003 

0.8 

0.9 

-0.13*  (0.01)  1 

8-9  Aug  89 

Palagic 

0.0061 

4.8 

Littoral 

0.0017 

4.7 

4.6 

-0.94*  (0.08)  { 

22-23  Aug  89 

Palagic 

0.0022 

2.0 

Uttoral 

0.0009 

2.3 

2.1 

-0.22*  (0.03)  1 

'  OMtarant  htttra  indieata  aignificafit  diflaranoM  at  itia  6-paroant  lawal  of  probaUBty  or  laaa 
(Ounean'a  muhipla  ranga  analyaia,  SAS  19861. 


were  greatest  in  early  July  and  early  August,  coinciding  with  the  occurrence 
of  relatively  large  mean  negative  horizontal  temperature  gradients  (Table  1). 
The  lowest  hourly  volumetric  flow  rates  were  observed  in  late  June  and  late 
July,  coinciding  with  minimal  mean  negative  horizontal  temperature  gradients. 
A  strong  linear  relationship  existed  between  hourly  volumetric  flow  rates  and 
mean  negative  temperature  gradients  (flow  =  4.38  ‘mean  negative  horizontal 
tenqierature  gradient  +  0.44;  r^  =  0.86;  p  <  O.OOS)  (SAS  1985),  suggesting 
that  differoices  in  water  temperature  between  littoral  and  pelagic  zones 
accounted  for  most  of  the  variation  in  houriy  volumetric  flow  rates. 


P  Exchange  Dynamics 

Hourly  TP  exchange  rates  to  both  the  littoral  and  pelagic  zones  w^-e  great¬ 
est  in  eariy  July  and  early  August  and  lowest  in  June  (Table  2).  TP  exchange 
r^es  to  both  zones  were  a  fimaion  of  both  the  volumetric  flow  rate  and  TP  in 
the  r^ion  of  flow.  Estimated  hourly  TP  exchange  rates  to  die  pelagic  zone 
were  greater  than  to  the  littoral  zone,  as  the  result  of  higher  TP  in  the  littoral 
bottom  wators  (Table  2).  Mean  summer  TP  in  the  region  of  flow  in  the  lit¬ 
toral  zone  (1 18.0  mg  m  ^  ±  15.7  SE)  was  over  2  times  greater  than  in  the 
region  of  flow  in  the  pelagic  zone  (44.5  mg  m  ’;  ±2.8  SE).  Net  TP  flux  to 
the  pdagic  zone,  at  rates  ranging  between  0.12  and  0.43  mg  m'^  hr',  occurred 
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Figure  7.  Vertical  and  longitudinal  variations  in  water  temperature  C’C)  measured  at 
0600  hr,  and  ^e  petition  of  the  dye  cloud  measured  between  0600  and  1 200  hr  on 
14  June,  28  June,  13  July,  27  July,  9  August,  and  23  August  1989.  Shaded  areas  repre¬ 
sent  position  of  dye  cloud 

during  all  dye-tracking  investigations  (Table  2).  Using  an  overall  mew  vol¬ 
umetric  hydraulic  flow  rate  of  1,770  m’  hr'  (determined  from  the  relationship 
between  hourly  volumetric  flow  rates  and  a  mean  negative  horizontal  temper¬ 
ature  gradiwt  of  0.41  **C)  and  mean  TP  in  the  region  of  flow,  a  mean  net  TP 
flux  of  0.22  mg  m'^  hr'  to  the  pelagic  zone  was  calculated  for  the  midsummer 
period. 
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1  Table  2 

1  Estimation  o 
1  Pelagic  Zone 
1  Tracking  Per 

■  ■  . . - . . . - . 

f  Hourly  TP  Exchange  Rates  Between  Littoral  and 

IS,  and  Net  TP  Flux  to  Pelagic  Zone  During  Dye- 
iods 

1 

Net  TP 

Houriy  TP 

Flux  to 

Diraetian  of 

TP 

Exchange  Rata 

Pelagic  Zona' 

Date 

Movamant 

mg  m"* 

mg  m'*  hr"' 

mg  m'*  hr’’ 

13-14  Jun  89 

Pelagic 

91.5 

0.22 

Uttoritf 

40.3 

0.10 

0.12 

27-28  Jun  89 

Pelagic 

124.8 

0.22 

1 

Littoral 

35.8 

0.06 

0.18  1 

1  12-13  Jul  89 

Pelagic 

134.6 

0.71 

1 

1 

Littoral 

52.8 

0.28 

0.43 

1  26-27  Jul  89 

Pelagic 

184.3 

0.27 

1 

Littoral 

43.8 

0.06 

0.21 

1  8-9  Aug  89 

Pelagic 

87.0 

0.67 

1 

Littoral 

52.2 

0.40 

0.27 

1  22-23  Aug  89 

Pelagic 

85.6 

0.30 

1 

Littoral 

42.0 

0.15 

0.15 

1  '  Oiffarartoa  batwaen  hourly  TP  oxchanga  to  tha  Httoral 

and  pelagic  zonae  d 

ividad  .  f  the 

_ 
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4  Discussion 


Rates  of  P  release  from  littoral  sediments  in  macrophyte  beds  of  Eau  Galle 
Reservoir,  measured  nr  the  laboratory;  were  linearly  i^ated  to  pH  under  both 
aerobic  and  anaerobic  conditions.  As  in  our  investigation,  other  studies  have 
shown  that  pH  can  have  a  regulating  effect  on  rates  of  P  release  from  sedi¬ 
ments,  particularly  under  aerobic  conditions  (Lee,  Sonzogni,  and  Spear  1977; 
Drake  and  Heaney  1987).  Enhanced  ligand  exchange  on  iron  hydroxide  parti¬ 
cles  at  high  pH  is  thought  to  be  an  important  mechanism  in  elevated  rates  of  P 
release  from  sediments  under  aerobic  conditions  (Drake  and  Heaney  1987). 

At  lowe^*  pH  values,  however.  P  release  from  littoral  sediments  appears  to  be 
minimal  under  aerobic  conditions.  This  pattern  may  be  attributed  to  the 
potentially  more  important  effects  of  an  anaerobic  environment  on  P  release  at 
lower  pH.  Estimated  rates  of  P  release  from  littoral  sediments  at  high  pH 
under  aerobic  conditions  were  in  agreement  with  diffusional  fluxes  estimated 
independently  from  sediment  peq>ers,  r^resent'mg  44  percent  of  the  mean 
reservoir-wide  internal  P  loading  rate  of  8  mg  m'^  day',  measured  for  the 
summer  by  mass  balance  in  Eau  Galle  Reservoir  (James,  Barko,  and  Taylor 
1991). 

TP  gradients  in  the  littoral  macrophyte  beds  were  prominent  dirou^out 
much  of  the  summer,  and  SRP  accounted  for  a  large  percentage  of  littoral  TP. 
To  our  knowledge,  no  investigations  have  shown  such  marked  concentration 
gradimts  vertically  in  the  littoral  zone  or  horizontally  between  littoral  and 
pdagic  zones.  These  observations,  in  combination  with  results  obtained  in 
our  laboratory  incubation  systems  and  sediment  peters,  strongly  suggest  that 
littoral  sediments  in  macrophyte  beds  provide  an  important  source  of  P  to  the 
resovoir.  Lower  P  concentrations  in  die  surface  watm  of  the  littoral  zone 
may  be  the  result  of  P  uptake  from  the  water  by  Ceratophyllum  demersum, 
since  the  foliage  of  this  species  is  most  dense  near  the  water  surfree  (James 
and  Barko,  personal  observation).  In  addition,  strong  thermal  stratification  in 
the  littoral  zone  during  the  day  (James  and  Barko  1991,  ui^ublished  data) 
may  restrict  upward  movement  of  sediment-derived  P  to  surface  waters  via 
mixing. 

Nighttime  convective  circulation  constituted  an  important  mechanism  for 
movement  of  littoral  P  to  the  pelagic  zone  during  the  summer.  During  the 
nigjit  and  morning,  cooler  bottom  water  from  the  littoral  zone  moved  as  an 
interflow  current  into  die  pelagic  zone  at  the  base  of  the  qiilimnion,  and  was 
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r^laced  by  a  return  current  of  wanner  surfoce  water  from  the  pelagic  zone. 
Tlie  hourly  volumetric  flow  rate  was  strongly  related  to  the  mean  negative 
horizontal  temperature  gradient  between  the  littoral  and  pelagic  zone,  indi¬ 
cating  that  convective  cooling  was  responsible  for  this  hydraulic  exchange 
patton.  Reasons  for  seasonal  variations  in  the  mean  negative  horizontal 
temperature  gradient  are  complex,  but  appear  to  be  related  to  the  rate  of  air 
temperature  cooling,  wind  effects,  relative  hunudity,  and  die  equilibrium 
tenqierature  (Stefan,  Horsch,  and  Barko  1989).  The  frequency  of  occurrrace 
of  negative  horizontal  temperature  gradients,  and  therefore  ni^ttime  convec¬ 
tive  circulation,  was  lowest  in  June,  but  increased  substantially  in  July  and 
August.  These  seasonal  variations  are  of  clear  importance  when  considering 
the  magnitude  of  littoral  P  flux  as  affected  by  nighttime  convective  circulation. 

We  have  focused  exclusively  on  nighttime  convective  circulation  as  a 
mechanism  transporting  littoral  P  in  an  underflow  to  die  pelagic  zone  of  Eau 
Galle  Reservoir.  However,  Monismith,  Imberger,  and  Morison  (1990)  have 
shown  recendy  for  another  small  reservoir  that  a  convective  overflow  pattern 
can  occur  during  daytime  heating  as  warmer  littoral  water  moves  as  a  surface 
flow  into  the  pelagic  zone.  Based  on  our  analysis  of  thermal  data,  these 
daytime  overflow  patterns  may  occur  in  Eau  Galle  Reservoir  as  well  (James 
and  Barko,  unpublished  data).  However,  since  surface  concentrations  of  P  in 
the  littoral  zone  were  typically  low,  and  often  similar  to  those  measured  in  the 
pelagic  qiilimnion  (Figure  4),  we  maintain  that  net  TP  flux  to  the  pelagic  zone 
of  Eau  Galle  Reservoir  during  the  day  is  probably  minimal. 

During  die  current  investigation,  it  was  assumed  that  hydraulic  and  P 
exchanges  determined  at  a  single  location  in  Eau  Galle  Reservoir  were  rqire- 
sentative  of  the  entire  littoral  zone.  However,  variations  in  the  slope  of  die 
littoral  zone  may  influence  the  rate  of  water  cooling,  because  of  differences  in 
volume  between  littoral  and  pelagic  zones.  Differences  in  slope  may  also 
affect  the  rate  of  bottom  water  movement  from  the  littoral  zone,  as  a  result  of 
differences  in  gravitational  acceleration.  The  effects  of  differing  macrophyte 
densities,  both  vertically  and  horizontally  in  the  water  column,  on  hydraulic 
exchange  rates  also  ne^  to  be  consider^. 

Prentki  et  al.  (1979)  estimated  a  net  P  flux  rate  of  3.2  mg  m'^  day’  for 
Lake  Wingra,  due  to  wind-induced  circulation.  However,  since  the  northwest 
bay  (region  of  study)  of  Eau  Galle  Reservoir  is  well  protected  from  northern 
and  western  winds  by  surrounding  hills,  and  is  densely  populated  by  C. 
demersum  (Filbin  and  Barko  1985;  James  and  Barko,  personal  observation), 
the  effects  of  wind  on  the  circulation  patterns  examined  in  this  reservoir  are 
consido'ed  to  be  relatively  minor.  Obviously,  this  may  not  ^ply  to  other 
systems  (e.g.,  Lake  Wingra)  with  greater  fetch,  surface  area,  or  exposure  to 
winds. 

Although  nighttime  convective  circulation  has  been  identified  as  an  impor¬ 
tant  transport  mechanism  for  littoral  P  in  Eau  Galle  Reservoir,  it  is  not  possi¬ 
ble  to  calculate  precisely  the  daily  net  TP  fluxes,  since  we  have  no  accurate 
measure  of  the  period  over  which  nighttime  convective  circulation  occurred. 
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In  earlier  investigations,  the  linear  movement  of  dye  in  this  reservoir  was  used 
as  a  measure  of  nighttime  convective  circulation  time  (James  and  Barko  1991); 
however,  this  information  was  not  obtained  in  the  present  study.  It  is  possi¬ 
ble,  however,  to  roughly  determine  the  period  of  nighttime  convective  circula¬ 
tion  as  the  time  during  which  negative  horizontal  temperature  gradients 
occurred.  This  period  averaged  9.6  hr  (±0.9  SE)  during  midsummer,  result¬ 
ing  in  an  estimated  mean  daily  net  TP  Oux  of  2. 1  mg  m'^  day  *  to  the  pelagic 
zone.  This  flux  rate  represents  about  60  percent  of  the  estimated  rate  of  P 
release  from  littoral  sediments  under  aerobic  conditions  (3.6  mg  m'^  day  ')  and 
26  percent  of  the  reservoir-wide  internal  TP  loading  rate  for  the  summer 
(James,  Barko,  and  Taylor  1991). 
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5  Conclusions  and 
Recommendations 


On  a  daily  basis,  the  water  in  macrophyte  beds  of  shallow  littoral  zones 
heats  and  cools  more  r!9)idly  than  adjacent  open  water,  partly  as  the  result  of 
differences  in  mixed  volume,  but  also  because  of  die  presoice  of  vegetation. 
Submersed  aquatic  macrophytes  contribute  significantly  to  the  development  of 
thermal  gradients,  which  promote  convective  hydraulic  circulation.  Effects  of 
convective  hydraulic  circulation  are  potentially  far-readiing,  since  this  mecha¬ 
nism  operates  daily  and  can  account  for  massive  transport  of  dissolved  constit¬ 
uents  between  littoral  and  open-water  regions  of  lakes  and  reservoirs.  In  Eau 
Galie  Reservoir,  as  reported  here,  convective  transport  from  aquatic  macro¬ 
phyte  beds  in  the  littoral  zone  accounted  for  about  26  percent  of  the  internal 
lo^  of  phosphorus  to  the  reservoir. 

For  other  aquatic  systems,  detailed  did  obsovations  of  both  water  temper¬ 
ature  and  periods  of  flow  wilt  be  necessary  to  estimate  convective  circulation 
time,  whi^  is  cr.tical  to  the  accurate  estimation  of  daily  horizontal  TP  fluxes. 
Investigations  of  hydraulic  aiKl  P  exchanges  during  convective  circulation  need 
to  be  conducted  under  a  variety  of  morphometric  and  climatic  conditions  to 
evaluate  for  other  aquatic  systems  the  significance  of  littoral  P  transport 
shown  to  be  important  in  Eau  Galie  Reservoir.  lu  addition,  it  is  important  to 
evaluate  the  effects  of  different  types  and  amounts  of  submersed  aquatic  mac¬ 
rophytes  on  convective  circulation.  This  latter  point  is  particularly  relevant, 
since  unlike  die  littoral  zone  of  Eau  Galie  Reservoir  (dominated  by  die  non- 
rooted  Ceratophyllum),  most  littoral  zones  are  dominated  by  rooted  aquatic 
macro{diytes. 

Information  on  rates  and  volumes  of  water  exchanged  in  littoral  macro¬ 
phyte  beJs  as  a  result  of  convection  and  other  processes  will  be  useful  in 
estimating  transport  and  dilution  rates  for  a  variety  of  solutes,  including  herbi¬ 
cides.  For  herbicides  specifically,  decisions  on  herbicide  type,  formulation, 
concentration,  and  time/season  of  application  will  be  direcdy  benefited.  This 
information  will  also  be  useful  in  assessing  the  influeices  of  aquatic  macro¬ 
phytes  on  eotro[diication. 
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